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We propose a method to obtain Floquet states — also known as light- induced states — and their 
quasi-energies from real-time wavefunctions without solving the Floquet equation. This is useful for 
the analysis of various phenomena in time-dependent quantum dynamics if the Hamiltonian is not 
strictly periodic, as in short laser pulses, for instance. There, the population of the Floquet states 
depends on the pulse form and is automatically contained in the real-time wavefunction but not in 
the standard Floquet approach. Several examples in the area of intense laser-atom interaction are 
exemplarily discussed: (i) the observation of even harmonics for an inversion- symmetric potential 
with a single bound state; (ii) the dependence of the population of Floquet states on (gauge) trans- 
formations and the emergence of an invariant, observable photoelectron spectrum; (iii) the driving 
of resonant transitions between dressed states, i.e., the dressing of dressed states, and (iv) spec- 
tral enhancements at channel closings due to the ponderomotive shift of above-threshold ionization 
peaks. 

PACS numbers: 32.80.Rm,02.70.Hm,32.80.Wr 



I. INTRODUCTION 

The time-dependent Schrodinger equation (TDSE) 
with a time-periodic Hamiltonian has solutions which can 
be expressed in a time-periodic basis. This basis is re- 
ferred to as the Floquet basis, and eigenstates in this 
basis are the Floquet states [IHJ]. Time-periodic poten- 
tials naturally arise when matter is exposed to laser fields. 
In this context, Floquet states are also known as "light- 
induced states" (LIS) [H, as they are the new states of the 
combined system "target -|- laser field". In fact, Floquet 
theory has been used to determine, e.g., very accurate 
ionization rates [gI, . Using so-called i?- matrix Floquet 
theory, the method has been extended to multi-electron 
systems [s^]. Strict periodicity of the Hamiltonian with 
the laser period implies physically that the laser pulse 
was always on and will be on forever. Then the prob- 
lem arises how the field- free system under study, e.g., an 
atom, gets into the laser field in the first place and how 
the field- free observables emerge. In fact, the population 
of the Floquet states depends on the laser pulse form. If 
the (up and down) ramping of the laser field is adiabatic 
and the laser frequency is non-resonant we expect the sys- 
tem to follow just a single Floquet state, namely the one 
which is adiabatically connected to the field-free initial 
state. However, for non-adiabatic ramping or resonant 
interactions a superposition of Floquet states is created. 
An example for non-adiabatic population of several Flo- 
quet states, leading to an apparent generation of even 
harmonics in inversion-symmetric potentials, is given in 
Sec. Uniof this paper. 

Instead of converting the TDSE into the time- 
independent Floquet equation [eq. ([13]) below] one may 
alternatively solve it directly in real-time. In the latter 
case there are no assumptions about periodicity or adia- 
batic ramping and, e.g., the effect of different laser pulse 
forms can be studied. However, the direct solution of the 
TDSE in real-time does not involve the Floquet basis 



so that information about LIS are not directly available. 
As many interesting phenomena such as the AC Stark ef- 
fect, Rabi oscillations, or stabilization against ionization 
[il, [l0| is most conveniently analyzed in terms of LIS, it 
is desirable to extract the "Floquet information" from 
the real-time wavefunction "on-the-fly" while propagat- 
ing (or by post-processing) it, without having to solve the 
Floquet equation as well. We present such a method to 
analyze non-perturbative, laser-driven quantum dynam- 
ics via the (time-resolved) Floquet information contained 
in the corresponding real-time wavefunction. 

The paper is organized as follows: in Sec. [Ill we re- 
view the basics of Floquet theory. In Sec. IIIII we briefly 
summarize the general derivation of harmonic generation 
selection rules before we present the (at first sight sur- 
prising) presence of peaks at even harmonics of the laser 
frequency in the case of an inversion-symmetric poten- 
tial with only one bound state. In Sec. |TVl we intro- 
duce our method to obtain the Floquet information from 
the real-time wavefunction, e.g., the populated states 
and their energies, and use them to explain the pres- 
ence of hyper- Raman lines at even harmonic frequencies. 
In Sec. |Vl we investigate how the population of Floquet 
states changes under (gauge) transformations while the 
Floquet energies and the observable photoelectron spec- 
tra remain invariant. In Sec. IVII time-resolved Floquet 
spectra of real-time wavefunctions in the so-called ve- 
locity gauge and in the Kramers-Henneberger frame-of- 
reference are compared. In Sec. IVIII the channel-closing 
phenomenon and related spectral enhancements are in- 
terpreted in terms of Floquet state-crossings. A conclu- 
sion is given in Sec. IVIIII In this work we restrict our- 
selves to spatially one-dimensional (ID) model Hamil- 
tonians. It is straightforward to extend the method to 
higher dimensions, as indicated in Appendix [Al Atomic 
units (a.u.) |e| = rrie = h = 47reo = 1 are used unless 
noted otherwise. 
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II. BASIC THEORY 

Consider a linearly polarized laser field E{t) of fre- 
quency uji in dipole approximation, polarized along the 
x-direction and interacting with an electron in some bind- 
ing potential V. The Hamiltonian in length gauge reads 



H{t) = Ho^ W{x, t), W{x, t) = xE{t) (1) 



with 
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A. Floquet theory 

For sufficiently long laser pulses 



(2) 



(3) 



holds to high accuracy, and thus also W{t + T) = W{t) 
so that 



H{t + T) = H{t). 



(4) 



The Floquet theorem states that in this case the 
TDSE 



i^*(x,t) = ^(t)*(x,t) 



(5) 



has solutions of the form 

^{x,t)=e-'''^{x,t), (6) 

^{x^t) being periodic itself, 

= ^(x,t + T). (7) 

6 is called the quasienergy or Floquet energy. The wave- 
functions ^{x^t) fulfill the Schrodinger equation 



n{t)^{x,t) = e^{x,t) 



(8) 



For a monochromatic laser field the interaction Hamilto- 
nian W{x^ t) can be written as 

W{x,t) = W'+(x)exp(icjit) + W'"(x)exp(-icJit), (12) 
leading to the time-independent Floquet equation 

(e + nfyjJi - Ho)(pn{x) (13) 

= W^{x)(pn^i{x) ^W-{x)ipn-l{x). 

The index n of the Floquet state is known as the "block 
index," which may be interpreted as the number of pho- 
tons involved in the process under study. Hence, the 
Floquet equation ([13]) couples any Floquet block n with 
its neighboring blocks nil via absorption or emission of 
a photon. 

In principle, (p!3|) is an infinite-dimensional set of dif- 
ferential equations. In practice, it is truncated so that 
^min < ^ ^ ^max- In obtaining the eigenvalue equation 
([13]), we assumed strict time-periodicity, which physically 
means that the laser pulse is always on. 

B. Non-Hermit ian Floquet Theory 

We are interested in systems which, in the field-free 
situation, possess besides bound states also a continuum. 
In the presence of a laser field such a system may ion- 
ize, i.e., the field- free stationary states are turned into 
field-dressed, quasi-stationary states. The simplest cases 
of only a few (field-free) bound states (allowing for res- 
onances) plus a continuum dressed by laser fields have 
been discussed in the literature since long ago (see, e.g., 
prj and [m for a review) . In an actual implementation of 
Floquet theory, the decay of quasi-stationary states needs 
to be taken into account when solving ([T3|) by applying 
Siegert boundary conditions for the outgoing waves 0, 
leading to complex Floquet energies 



e = Re e — i- 



(14) 



with 



where F is the ionization rate. The difference between 
Ree and the field-free e^^^ is the AC Stark shift. 



C. Finite-grid, finite-pulse TDSE solution 



n{t) = H{t) - i-. 



(9) 



If e is an eigenvalue and ^{x^t) the corresponding eigen- 
state, also 

e' = e + ma;i, <I>'(x, t) = e^^^^^^(x, t), mGZ (10) 

are solutions of ([HJ. Owing to the time periodicity of 
$(x,t) we can expand 



\noj-it 



(11) 



We solve 



i^*#(x,t) = ff(t)*#(x,t) 



(15) 



§ < X < -| for times 



on a numerical grid of size L, 

< t < tsim with tsim the total simulation time. The 
binding potential V{x) is centered at x = 0. In all cases 
discussed in this work we start from the field-free ground 
state on the grid ^^(x,0) = ^^^(x). Probability den- 
sity approaching the grid boundary is absorbed by an 
imaginary potential. 

Our aim in the following will be to analyze ^^{x^ t) in 
terms of Floquet energies and states. 
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III. HARMONIC GENERATION 

In the first example we apply our method to investigate 
the origin of apparently even harmonics in an inversion- 
symmetric potential with only one bound state. 

There are many ways to derive selection rules for har- 
monic generation (HG). Most elegant, rigorous, and ap- 
propriate for our purpose is the approach employing dy- 
namical symmetries [H, Consider the stationary 
Schrodinger equation 
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(16) 



with Hq given by ([2]). If the potential V is inversion- 
symmetric, V{x) = V{—x), the Hamiltonian Hq is in- 
variant under spatial inversion as well. 



Ppfix) = f{-x), P2 ^ 1^ p- 



(17) 



(18) 



SO that for non-degenerate energies £ the eigenstate ^{x) 
is also an eigenstate of the spatial-inversion operator Pp. 
Because of = 1 the eigenvalues can only be ±1 (par- 
ity): 



Pp^(x) = ±^(x). 



(19) 



The full Hamiltonian ([T]) and the Floquet-Hamiltonian 
(|9j) are not invariant under spatial inversion but under the 
dynamical symmetry operation "spatial inversion com- 
bined with a translation in time by half a period" , 



[^W,^pt] = [H(t),Ppt]=0, 



Pptf{x, t) = /(-X, t + tt/cji). 



p2 



1, 



(20) 



(21) 



pt 



PpPt 



PfPp, 



Ptf{x,t) = f{x,t+7r/0Ji). (22) 



For non-degenerate e 

PpMx,t) = ±^x,t). 
Because of (pT]) we observe that 

PpMx,t) = ^(-l)"e-'"'^^*4^„(a;), 



and with 



follows that 

PpLPn{x) = ±{-l)"(pn{x), 



(23) 
(24) 

(25) 



i.e., the (pn{x) have an alternating parity with respect to 
the Floquet block index n. 

Numerically, the HG spectrum ~ a;^|(i(a;)p is calcu- 
lated via the Fourier-transformed dipole moment 



d{uj) 



+ L 

J dt J dx^^^{x,t)x^^ 



(x, t) e 



iujt 



(26) 
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FIG. 1: (color online). Logarithmically scaled HG strength 
uj^\d{uj)\'^ vs harmonic order and excursion amplitude a — 
A/(jJi {(jJi = 1, vector potential A{t) ramped up and down over 
4 cycles and held constant with amplitude A for 30 cycles). 
The numerical fast-Fourier transform was performed over the 
pulse duration, i.e., tsim = 38 cycles, using a Hanning window. 



Assuming that the numerically determined, exact wave- 
function on the grid is well described by just a single 
Floquet state, using (j6j), (pTj) , and (p!4|) yields 



d{uj) ip*^{x)xipn{x) dx 



(27) 



/' 

Jo 



^t{i[u;-u;i(n-m)]-r} 



The spatial integral is non- vanishing only if (fn has the 
opposite parity of cpm^ i-e.. 



n — m = 2/c + 1, 



kez. 



(28) 



The temporal integral thus leads to peaks centered at 
frequencies 



(29) 



with widths determined by tgim (frequency-time uncer- 
tainty) and F (decay). The selection rule (|29|) is the 
well-known result that an inversion-symmetric target in 
a linearly polarized laser field generates odd harmonics 
only. Note that the above derivation also holds for multi- 
electron targets because the electron-electron interaction 
is also invariant under the symmetry operations Pp and 

Ppt- 



A. Hyper-Raman lines at even harmonics of the 
laser frequency 

It is known that HG peaks at positions different from 
odd multiples of the fundamental laser frequency Ui are 
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to be expected for an inversion- symmetric potential if at 
least two Floquet states of opposite parity are populated 
Physically, the superposition of two Floquet 
states may amount to, e.g., the absorption of n photons of 
energy uji but emission of one photon of energy nuji — Ae, 
with Ae being the energy difference between initial and 
final state. This should lead to hyper-Raman lines in the 
spectra which, however, are typically weak [ill, U^- Nev- 
ertheless, if observable, they appear at even harmonics of 
the laser frequency in the case of degeneracy, Ae = 0. 
We consider an electron in the Poschl- Teller potential 



V{x) 



1 



cosh X 



(30) 



and subject to a laser field. The potential (|3Q|) supports 
only a single bound state ^o(^) of energy So = —0.5. 
Hence, superpositions of field-free bound states are ruled- 
out. As a consequence, perturbation theory in the exter- 
nal field can certainly not predict hyper-Raman lines or 
even harmonics. However, Fig. [T] shows the logarithmi- 
cally scaled HG strength as obtained from the 
numerical solution of the TDSE. The HG strength is plot- 
ted vs harmonic order uj/uoi and the amplitude a of the 



excursion 



a{t) 



A{t) dt 



(31) 



with A{t) the vector potential of the laser field. The 
electric field is given by E{t) = —dtA{t). Given the 
vector potential amplitude A, the excursion amplitude 
is a = A/cJi, the field amplitude E = Auji. The laser 
pulse parameters are specified in the figure caption. One 
sees that for sufficiently strong excursion amplitude a 
peaks at even harmonics of the laser frequency appear 
too. Picking an even harmonic at a > 15 (e.g., the 6th) 
and tracing it back to low a reveals that the peak splits 
and rapidly drops in magnitude (e.g., around a 2 for 
the 6th harmonic). In the next Section we will use our 
real-time Floquet method to show that the appearance 
of even harmonics is due to the population of several LIS 
that become quasi-degenerate as a increases. 



B. Superposition of Floquet states 

In the case of a non-adiabatic transfer of the field-free 
state to field-dressed states one has to allow for a super- 
position of Floquet states in order to represent the exact, 
numerically determined wave function on the grid, 

(32) 

Here we assume that the expansion coefficients are in- 
cluded in $^(x, t) and (p^ri{x). For continuous quasiener- 
gies the sum over /3 should be replaced by an integral over 
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FIG. 2: (color online). Logarithmic plot of = |Q+|^ + |Q_ |^ 
vs energy 8 and excursion amplitude a = A/uji, showing the 
quasi-energies of the (populated) field-dressed states. The 
laser frequency was cji = 4. The pulse shape was trapezoidal 
(4,1200,4) in the vector potential of amplitude A. For each a 
the maximum in R was renormalized to unity. 



e. The Fourier- transformed dipole will be 



d{uj) ^ / (p*^{x)xip^n{x)dx 



(3-ynm 



(33) 



Jo 



^^{i[u;-u;i(n-m)-(Ree^-Ree^)]-(r;3+r^)/2} 



Again, in order for the spatial integral to not vanish 
the parity of cpf^n and cp^rn must be different. However, 
now this can be the case not only for n — m = 2/c + 1, but 
also for n — m = 2k if the parity of, e.g., (piso is opposite to 
the one of (p^Q. Hence, one expects the above-mentioned 
hyper-Raman peaks at 



kuji + Ae, 



(34) 



where Ae = Ree/3 — Ree^ is the difference between 
the real parts of the two Floquet quasienergies involved. 
Thus, in order to observe even harmonics at exactly 
u = 2kuoi a degeneracy Ree/5 = Ree^ is required. Such 
a degeneracy between the (field-dressed) initial state and 
another one of opposite parity is also likely to populate 
the latter one. 



IV. FLOQUET STATE ANALYSIS OF 
REAL-TIME WAVEFUNCTIONS 

The extraction of Floquet information contained in the 
real-time wavefunction is useful to analyze any feature of 
interest in HG spectra. We start with the determination 
of the (real part of the) quasi-energy of the populated 
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FIG. 3: Field- dressed ground state wavefunction (pon for a — 
4. (a) Floquet block n = 0, (b) n = -1. 

Floquet states. Once these energies are known the corre- 
sponding Floquet states can be obtained. The method is 
similar to the one proposed in pisj for field-free dynamics. 

The numerical solution of the time-dependent 
Schrodinger equation in real time yields Upon 
multiplication of (|32]) by an even or odd test function 
q±{x)^ spatial integration, and Fourier transformation 
from the time to the energy domain, 

t — 

(35) 

< ti < t2 < tsim one can extract from the peak posi- 
tions in \Q±{E)\'^ the real part of the Floquet energies 

Re Ei3n = Re 6/3 + nui (36) 

belonging to even or odd Floquet states (pisn^ respec- 
tively. The even test function is, e.g., simply unity for all 
— ^ < X < ^, the odd test function may be chosen 1 for 
X > and — 1 for x < 0. The purpose of these test func- 
tions is to extract the even and odd parity Floquet states 
separately. Of course, only the energies of the populated 
(and thus relevant) Floquet states cpf^n are obtained in 
this way. 

The imaginary part ofe^ contributes to the width 
of the peaks in |Q±(f)p. However, in our finite-time. 



finite-grid TDSE simulations the width of the peaks in 
\Q±{£)\'^ also depend on the integration time t2 — ti and 
the grid-size because of the absorbing grid boundaries. 
Only for a flat-top laser pulse and a very long simulation 
time a stationary absorption rate at the grid boundaries 
would be estabhshed, and could be determined from 
the peak-width. This, however, is exactly the regime 
where the standard Floquet approach should be applied. 
We focus here on aspects of our method complementary 
to the conventional Floquet method, in particular its ap- 
plicability to finite pulses and time-resolved studies. 

If we multiply the wavefunction (|32]) by exp(itf ) (with 
£ real) and integrate over time, mainly the Floquet state 
(fs for which the phase is stationary, i.e., £ = Re e/3 + ncJi 
"survives," 

ips{x) - [ \''^^{x,t)dt. (37) 

The integration time t2 — h has to be sufficiently long 
in order to cover many temporal oscillations of the wave- 
function. 

Starting from the ground state in the potential ([3Q|) . we 
solved the TDSE for a high-frequency laser field of vec- 
tor potential A{t) = —A{t) sinujit for cji = 4 and A{t) a 
trapezoidal pulse shape with linear up- and down-ramps 
over 4 cycles and 1200 cycles constant amplitude A (de- 
noted in the form (4,1200,4) in the following). Figure [2] 
shows 

R = \Q+f + \Q.f (38) 

(with the time-integral in (|35]) performed over the en- 
tire pulse) as a contour plot vs the excursion amplitude 
a = A/uji and energy £ for an energy interval within 
the zeroth Floquet block n = 0. Plotting and 
IQ-P individually allows to distinguish the parity of the 
states (labeled 'even' or 'odd' in Fig. [2j). For a ^ 
only the field- free state at £ = —0.5 remains. How- 
ever, with increasing excursion amplitude a light-induced 
quasi-bound states emerge, which are populated due to 
the finite rise-time of the laser field. From the popu- 
lations (see color-coding) one infers that around a = 6 
besides the field-dressed ground state the second excited 
field-dressed state is more populated than the first ex- 
cited. For increasing a the field-dressed ground state 
and the field-dressed first excited state become almost 
degenerate so that Ae ^ in explaining the peaks 
at even harmonics of the laser frequency due to hyper- 
Raman scattering. 

Using (|37|) we extracted field-dressed states. Fig- 
ure [3] shows the field-dressed ground state for the Floquet 
blocks n = (a) and n = — 1 (b) for a = 4. The integra- 
tion time was again the pulse duration. Equation ([37|) 
in general yields a complex wavefunction (ps = (ps^^^ - 
The plots in Fig. [3] show the real wavefunction (p)^. It is 
seen that the parity indeed changes as one decreases n 
by one. For n = and a = the ground state must be 
even. Hence, for n = —1 it is odd, in accordance with 

(ESI). 
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FIG. 4: (color online). R vs energy 8 and second-laser 
frequency Cu for first-laser excursion a = 2.5 and second- laser 

field strength amplitude E = 0.01. 



V. TRANSFORMATIONS 

We consider transformations G{t) which are periodic 
in time and reduce to unity as the laser field goes to zero, 

G{t + T) = G{t), G{t)\a,E,A=o = 1. (39) 
Now, since each Floquet state fulfills (|8]), 

G{t)mG-\t)Gmm = n'm'^it)) = e^\^'^{t)) 

(40) 

where n\t) = G{t)n{t)G-^ (t) is the transformed 
Floquet-Hamiltonian and |<^>^(t)) = G{t)\^i3{t)) the 
transformed Floquet state. The quasi-energy is not 
affected by the transformation, and \^^^{t)) is also peri- 
odic because of (|39|) . so that with (pTj) 



It is known that if the laser frequency is tuned around 
resonances field-dressed states originating from different 
Floquet blocks (and corresponding to the coupled field- 
free states) display avoided crossings. These crossings 
have been shown to be related to localization, and to 
chaos in the corresponding classical system [19|. The 
separation of the two dressed states involved corresponds 
to the Rabi frequency and is proportional to the field 
strength of the driving laser. We will now show that the 
same is observed for transitions between already dressed 
states, i.e., we use the laser of frequency uji to dress 
the system and a second, weaker laser of frequency uj 
to induce transitions between dressed states. The second 
laser will dress the already dressed system [20|, and the 
"dressed^" states (or two color-dressed states) should dis- 
play avoided crossings as the frequency Cj is tuned around 
the energy gap of two dressed states. 

From Fig. [2] one infers that for an excursion ampli- 
tude, d = 2.5 the energy difference between the field- 
dressed ground state and the field-dressed first excited 
state is Reei — Reeo — 0.155. Hence, we tune the fre- 
quency UJ of the second laser around this energy differ- 
ence. The pulse envelope was the same for both lasers, 
and the electric field amplitude of the second laser was 

E = 0.01 = Auj = auj'^ for all uj. Figure |4] shows results 
for the Floquet energy spectrum R vs energy and uj for 
d = 2.5. If the two frequencies uji and uj are incommen- 
surate the Hamiltonian is not periodic at all. However, 
our approach does not require periodicity, and we expect 
a Floquet analysis to be meaningful as long as the two- 
color Hamiltonian is approximately periodic, namely in 
T = 27t/uj because uji ^ Co. In fact, the avoided cross- 
ings of Re eo with Re ei — uj and of Re eo + with Re ei 
around uj = 0.155 are clearly visible in Fig. HJ 



where G{t) = T.^e-^'^^'^Gm, and thus 



(42) 



We now specialize on transformations G that commute 
with the dynamical symmetry operation Ppt, 



[G{t),Ppt] = 0. 



(43) 



Examples are gauge transformations, e.g., for the trans- 
formation from velocity gauge, where 



1 



W{t) = pA{t) 

to the length gauge one has 

GlgW = exp[ia;A(t)]. 



(44) 



(45) 



Another example is the Pauli-Fierz or Kramers- 
Henneberger (KH) transformation, which is not a gauge 
transformation (although one frequently finds the term 
"KH gauge" in the literature). If we start from the veloc- 
ity gauge interaction dH]) the KH transformation reads 



GKuit) = exp 



f A\t')dt' ^\a(t)p 



(46) 



This amounts to a translation in position space by the 
free electron excursion a{t) (|3T]) and a purely time- 
dependent contact transformation. The KH Floquet- 
Hamiltonian is 

n'{t) = UKRit) = -p" + V[x + a{t)] - i-. (47) 

As a consequence of 

P^tWAt)) = G{t)P^t\^p{t)) = ±\<i>'M)) (48) 
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with the eigenvalue ±1 the same as for Ppt|$^(t)) = 



One also finds Gm = (— 1)^ 



p n 

1 pKjTri 



Prj and 



±{—lY\(p'^^)^ i.e., the transformed (primed) 



states have the same symmetry as the original states. 

Figure [5] shows the KH and the velocity gauge proba- 
bility density for the excursion amplitude d = 10. The 
target energy was 6 = —0.08 where in Fig. [2] the al- 
most degenerate ground and first excited state energies 
for a = 10 are. The KH probability density fits to the 
KH potential 



VKuix) 



1 



V[x -\- a sinr] dr, 



(49) 



shown in the lower panel. The actual calculation was 
performed for cji = 4 and a trapezoidal (10,1180,10)- 
pulse. The target energy £ in ([37|) is scanned through 
the energy region of interest, and the Floquet energy is 
hit when the value of the integral is maximum. If one uses 
the same integration time for different £ the integral 



t 



\ips(x)\^Ax 



(50) 



is a relative measure for the population of the respective 
Floquet state in the actual pulse. 

The Floquet energies are invariant under the transfor- 
mations Git) while both the Floquet states \^^r^ and 
their populations are not. In particular, in the high- 
frequency limit one expects that only the eigenstates in 
the KH potential (|49]) matter [9]. These states corre- 
spond to the Floquet energies in the Floquet block n = 0. 
Hence, the energy spectrum in the KH frame is expected 
to be much more localized around n = than in velocity 
gauge. This is confirmed by Fig. [6l Instead of using the 
even or odd test functions in ([35]) and spatial integration 
we analyzed the wavefunction ^^(x,t) at Xtest = 2, i.e., 
we calculated 



E 

/3n 



-it(e,+„c.,-£)d^ ^/J„(a;test). (51) 



This avoids the transformation of the entire wavefunc- 
tion to the KH frame and yields similar results as long 
as one chooses Xtest in a region where the wavefunction 
is sizable and both odd and even parity wavefunctions 
contribute (for Xtest = only contributions from even 
Floquet states would be visible). Figure [6] confirms that 
for transformations of the type (|39|) the populations of 
Floquet states in different frames (or gauges) are differ- 
ent while the Floquet energies are the same. The latter, 
dressed levels could be probed with a second laser [21^ . 
Of course, any gauge- or frame-dependence should vanish 
when field-free observables, such as photoelectron spectra 
are considered. 



VI. PHOTOELECTRON SPECTRA 
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FIG. 5: (a) KH and the velocity gauge probability density for 
the excursion amplitude a = 10 and target energy £ — —0.08. 
(b) Corresponding KH potential. 



Teller potential have energies f > 0. With laser field all 
continuum and bound states are contained in each Flo- 
quet block so that overlaps of dressed bound states from 
one block with continua from other blocks with lower 
n are possible. However, we expect the dressed bound 
states of the n = block to dominate since they are the 
main ones being populated during the switching-on of the 
laser. Let us first discuss the case where uj\ > minRee/3, 
i.e., a single photon is sufficient for ionization. Then 
the dressed bound state in Floquet block n with energy 
Re 6/3 + na;i overlaps with continuum states of all the Flo- 
quet blocks m < n. In particular, Ree/3 + nui overlaps 
with the continuum state of energy Cp of the zeroth Flo- 
quet block, where p indicates the asymptotic momentum 
of this continuum state. 

We will now turn to the question of how the manifold 
of mixtures of bound and continuum Floquet states con- 
verts to an observable photoelectron spectrum when the 
pulse is switched off. Figure shows a time-resolved Flo- 
quet spectrum in velocity gauge for a A^cyc = 100-cycle 
sin^ 



"^-pulse 



sinujit 



(52) 



Without laser field the continuum states of the Poschl- for < t < iVcyc27r/a;i and zero otherwise. The other 
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FIG. 6: Floquet spectra for a — 10, cji = 4, and a 
(10,1180,10)-pulse in (a) velocity gauge (with the A^{t)/2- 
term transformed away) and (b) in the KH frame. In the 
KH frame the n = 0-Floquet block dominates while in ve- 
locity gauge the population is broadly distributed over many 
Floquet blocks. 
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FIG. 7: (color online). Time-resolved Floquet spectra for a 
100-cycle sin^-pulse of amplitude a = A/ui = 10, uoi =4, 
attest = 2 (i.e., "inside" the potential), and a time- window of 
width tw = t2 — ti = 50. The vertical line indicates the end of 
the pulse. Panel (b) is a close-up of the energy region around 
eQ^^* = —0.5 in (a). The calculation was performed in velocity 
gauge (with the A^(t)/2-term transformed away). 



pulse parameters are given in the figure caption, and 
^test = 2 (i.e., "inside" the potential) and a time- window 
of width = t2 — ti = 50 were chosen for (|5T]) . The 
time on the horizontal axis is ti so that the spectrum for 
times ti > lOOT = 157.1 (indicated by the vertical black 
line) shows field- free states, i.e.. 



(53) 



/3 "^^1 



Figure [7^ shows that while the pulse is on the population 
is distributed over many Floquet blocks. As the pulse is 
switched off, all the Floquet populations for n 7^ dis- 
appear, and only the ground state population inside the 
potential with energy eQ^^ remains. This is because we 
analyzed the spectrum at the position Xtest = 2. Contri- 
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FIG. 8: (color online). Same as in Fig.[71but for xtest = 471.3. 
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butions to the wavefunction corresponding to electrons 
in the continuum, travehng with an asymptotic momen- 
tum decay at Xtest = 2. Figure [Tb shows a close-up 
of the region around eQ^\ With increasing amplitude of 
the laser pulse the dominant Floquet population shifts 
adiabatically from the field- free value eg^^ = —0.5 to the 

ground state energy of the KH potential e[)™^ -0.09 
(see Fig. [2] for d = 10) and back. Note that although 
the calculation was performed in velocity gauge the KH 
ground state energy is relevant here because the Floquet 
quasi-energies are frame- and gauge-independent. 

Figure [8] shows the same analysis for xtest = 471.3, i.e., 
"far away" from the atom so that it takes some time until 
probability density arrives there, namely around t = 100. 
It is interesting to observe that in velocity gauge this "ar- 
rival time" during the pulse is independent of the energy. 
As the laser pulse is switched off at t = 157.1 many Flo- 
quet channels close. However, because electrons are still 
on their way from the atom to the "virtual detector" 
at Xtest = 471.3 we are able to "measure" the field- free 
photoelectron spectrum of the electrons emitted in that 
direction. The time these free electrons need to pass the 
virtual detector decreases with increasing energy, as is 
seen in Fig.[8]where the width of the traces for t > 157.1 
decrease with increasing energy. The five traces visible 
are separated by uji and correspond to above-threshold 
ionization (ATI) peaks (see, e.g., the review (22j or (23|). 
They are quite broad in energy because of the change 
of the ionization potential (from field-free value to KH 
value and back). Their figure-eight shape in the contour 
plot of Fig. [8] is a peculiarity of the sin^-pulse shape. 

Figure [9] shows the corresponding result obtained in 
the KH frame. We see that in the KH frame only those 
states are populated in the laser field which actually con- 
tribute to the final field-free spectrum. This is because 
the KH potential at Xtest = 471.3 is almost identical to 
the field-free one so that outgoing electrons are not af- 
fected anymore by the oscillating KH binding potential. 
It is also seen in Fig. [9] that the most energetic electrons 
arrive earlier at Xtest, unlike the velocity gauge-result in 
Fig. El 



VII. CHANNEL-CLOSINGS 

So far we studied mainly high-frequency phenomena 
where the Floquet blocks are well separated on the 
atomic energy scale because the laser frequency exceeds 
the ground state ionization potential. However, there are 
plenty of interesting, non-perturbative phenomena occur- 
ring at low frequencies where the ponderomotive energy 
Up = E'^/Aujf can be large at nowadays available laser 
intensities E'^. Examples are tunneling ionization and 
high-order ATI due to rescattering of electrons [23|, [l^ . 
In this Section we choose the so-called "channel-closing" 
(see [25[ and references therein) as a low-frequency phe- 
nomenon to illustrate our method. 




Time (a.u) 

FIG. 9: (color online). Same as Fig. [8] but in the KH frame. 



The TDSE was solved for a trapezoidal pulse of fre- 
quency uji = 0.08. On the energy scale of the ioniza- 
tion potential the Floquet blocks are packed much closer 
in this case, meaning that many photons are necessary 
for ionization. In Fig. [10] we plot the Floquet energy 
spectrum in a certain range of excursion amplitude 
d = E/uOi and energy £ around the field- free contin- 
uum threshold (other relevant parameters given in the 
figure caption). The calculation was performed in ve- 
locity gauge using again the potential (|3Q|) . There is a 
clear down-shift of all the populated Floquet levels with 
increasing laser amplitude. This AC Stark shift is also 
referred to as the "ponderomotive shift" because the ef- 
fective ionization potential is increased by U^. In fact, 
the energy in the photoelectron spectrum is given by 

f = y =nnc^i-(|^o| + t/p), (54) 

(provided the AC Stark shift of the initial state is negligi- 
ble, which for atomic ground states at long wavelengths 
often is the case). £q is the initial electron energy and n 
is the number of photons absorbed. In order to reach the 
continuum at all n > (|£^o| + Up)/Jiwi photons have to 
be absorbed. As the intensity, and thus /7p, is increased, 
more and more photons are needed for ionization. When 
n photons are no longer sufficient but n + 1 photons are 
needed the n-photon channel closing occurs. In the plot 
shown in Fig. [TO] a channel closing manifests itself as a 
crossing of a Floquet quasi-energy and the continuum 
threshold. Now, the interesting feature in Fig. [10] is the 
zero-energy LIS. Such LIS were also observed in Ref. p6j . 
where their connection with experimentally observed en- 
hancements in the photoelectron spectra at high energies 
p7j was established. The parity of both states involved 
in the crossing in Fig. [10] is even, and it is known that 
depending on the parity of the states, channel closings 
affect the photoelectron spectrum differently [26|, [28| . 

In our model, for the first even channel closing eight 
photons are needed. According ([54]) it is expected at 
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FIG. 10: (color online). Logarithmic plot of = |Q+|^ + 
IQ-I^ vs energy 8 and excursion amplitude a, showing the 
(populated) field-dressed states for cji = 0.08 and a trape- 
zoidal (4,40,4)-pulse. 



a = 9.354, which indeed is close to where the crossing 
is observed in Fig. [TOl The small discrepancy is because 
of the AC Stark shift of the initial state, neglected in 
([5^ . One would expect that channel closings only affect 
low-energy electrons because the kinetic energy of the 
electrons whose channel is about to close is low. Hence, 
as the intensity is increased the yield of ATI peaks at 
energies, say, > 5[/p should increase monotonously as 
well. However, near even photon channel closing there 
is a marked increase in the photoelectron yield at high 
energies [25|, IH, |28| . Instead, when in odd photon chan- 
nel closings the odd-parity LIS crosses the zero-energy 
LIS, such enhancements are absent. The first odd pho- 
ton channel closing occurs around a = 11.55, the next 
even photon channel closing occurs around a = 13.55. 
The photoelectron spectra obtained using our Floquet 
method confirm the presence and absence of enhance- 
ments at even and odd channel closings, respectively, as 
shown in FigfTTl 



VIII. CONCLUSIONS 

We described a method for obtaining Floquet infor- 
mation from real-time wavefunctions. In this approach, 
it is not necessary to assume strict periodicity. In fact, 
it is possible to follow the time-resolved Floquet quasi- 
energies as they shift during a laser pulse. Moreover, 
the populations of the Floquet states can be determined 
so that especially cases where superpositions of Floquet 
states play a role can be identified. The usefulness of 
the method was illustrated by several examples employ- 
ing the one-dimensional Poschl- Teller potential with only 
a single field- free bound state. In particular, the ori- 
gin of peaks at even harmonics of the laser frequency 
in an inversion-symmetric potential, avoided crossings 
of dressed already field-dressed states induced by a sec- 
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FIG. 11: (color online). Photoelectron spectra around 5?7p, 
parameters as in Fig. 1101 (a) Non-monotonic behavior of the 
yield (open squares a — 13.0, solid squares 13.3, circles 13.55, 
triangles 13.8). (b) Same for an odd channel closing, showing 
a monotonic behavior of the yield with increasing intensity 
(triangles a = 11.8, circles 11.55, open squares a = 11.3). 



ond laser, the properties of Floquet states under time- 
periodic transformations, the emergence of invariant, ob- 
servable photoelectron spectra after the laser pulse, and 
photoelectron enhancements at channel closings were dis- 
cussed. The method is straightforwardly extendable to 
three dimensions. We think the method is most useful 
for researchers running codes to solve the time-dependent 
Schrodinger equation in real time. By saving the wave- 
function at selected spatial positions as a function of time 
during the interaction with the laser field the analysis in 
terms of light-induced states can be easily performed a 
posteriori. The application of the method to correlated 
multi-electron systems may be very fruitful, as the un- 
derstanding of field-dressed, multiply-excited or autoion- 
izing states is still poor. 
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Appendix A: Extension to three dimensions 

The method of Floquet analysis described in this work 
is easily extendable to higher-dimensional systems. For, 
e.g., hydrogenic systems in 3 dimensions (3D) one could 
follow the evolution of Floquet states with different or- 
bital angular momentum quantum numbers Instead of 
(|6|) we have 

^{r,e,cP,t) = e-'''<^>{r,e,cP,t) (Al) 
with [compare to (pT]) ] 

$(r, e, cP,t) = J2 V>n{r, e, 0)e-'""i*. (A2) 

n 

The operator Ppt [compare to ([2T]) ] acts according 

Pptf{r,t) = f{-r,t+—), (A3) 
and (|24|) becomes 
Ppt^{r, e, (/), t) = ^(-1)" exp(-inc^it)Pp(^n(r, 0, (/>). 

n 

(A4) 



If we expand the (fni'f^-, ^i^) in spherical harmonics. 



^n{r,0,^)=Rni{r)Yi^{0,^), (A5) 



we find, using 



PpYira{e, (/.) = Fz^(7r - ^, TT + 0) = {-l)'Yira{9, <^), (A6) 



that 



Pp^r, e, 4>) = (-i)"+Vn(r, e, 4>), (A?) 



the analogue of (|25]) . Note that n is the Floquet block 
index here, not the principal quantum number. After 
these considerations it is straightforward to extend the 
Floquet analysis of real-time wavefunctions described in 
this work to 3D. 
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